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ABSTRACT: In the present study, we have accounted for the accelerated expansion of the universe in the 
framework of the generalized Brans-Dicke theory. It has been shown very clearly that the transition of the universe 
from a phase of deceleration to that of acceleration can be explained by the generalized Brans-Dicke theory, without 
having to consider any contribution from the entity called dark energy. The present model shows a signature flip of 
the deceleration parameter, from a positive to a negative value, occurring in the recent past of the matter 
dominated era. We have determined the functional form of the Brans-Dicke parameter, expressed in terms of the 
scalar field. This functional form causes the effective pressure to be sufficiently negative to drive the accelerated 
expansion of the universe. It has been found in the present study that the gravitational constant increases with time. 
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1. INTODUCTION 

It has been established beyond doubt, by recent 
astrophysical observations, that there was a smooth 
transition of the universe from a decelerated phase 
to an accelerated phase of expansion [1, 2]. The 
reason behind this accelerated expansion has been a 
matter of intense research for many years. It is 
known to us that normal matter interacts through 
the gravitational force and it has a positive definite 
density and pressure. Therefore, there must be some 
other kind of matter responsible for the accelerated 
expansion, which makes the effective pressure 
sufficiently negative, causing a repulsive effect. This 
new entity is generally referred to as dark energy. In 
the past few years a lot of research has been 
undertaken to analyse the nature of this new form of 
energy, responsible for accelerated expansion. A 
number of models have been proposed so far to 
account for the acceleration of the universe. The true 
composition and behaviour of this new entity is yet 
to be determined. 

The most obvious choice to represent dark energy 
has been the cosmological constant, denoted by the 
symbol A (lambda) [3]. But there are serious 
problems regarding the value of this constant which 
can’t be accurately determined by the ACDM model. 
The value of A, obtained from recent observations for 
an accelerating universe, is not consistent with the 
value in Planck scale or Electro-weak scale [4]. By 
constructing a model, where A is a function of time, 
this problem can be partly solved. Several such 
models have been formulated but they have their 
own shortcomings [5, 6]. 

Scalar field models have been proposed as an 
alternative to the dynamical A models. Here, the 
equation of state of dark energy is a function of time. 
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Quintessence models were proposed as one of the 
most important types among the scalar field models. 
These models have a potential function V(@), 
enabling the pressure sector, which is a function of 
the scalar field, to have a sufficiently large negative 
value, leading to the predictions consistent with 
observations regarding the accelerated expansion of 
the universe [7, 8]. A major limitation of these 
quintessence models is that most of the quintessence 
potentials are arbitrarily chosen, without having any 
theoretical basis behind their formulation. Quite 
naturally, a large number of models, based on scalar 
fields, such as the tachyon, k-essence, holographic 
dark energy models have come into existence with 
their strengths and limitations [9, 10, 11]. 


The cold dark matter and dark energy are normally 
allowed to evolve independently in most of the scalar 
field models. There are some studies where one of 
them decays into the other through some interaction 
between them [12]. A non-minimal coupling of the 
scalar field with the dark matter sector, through an 
interference term in the action, has helped us to 
explain the accelerated expansion of the universe. 
These are known as chameleon fields and they can 
effectively serve as dark energy, causing the cosmic 
acceleration [13, 14]. 


In the framework of Brans-Dicke (BD) theory, non- 
minimal coupling between the scalar field and 
geometry, play a very important role in explaining 
the observed acceleration. A _ self-interaction 
potential V(@), which is a function of the BD scalar 
field @, has often been incorporated in the 
framework of Brans-Dicke theory to explain the 
cosmic acceleration [15]. In some cosmological 
models, a quintessence scalar field is introduced in 
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the BD theory. It can account for the late time 
acceleration for a wide range of potentials [16]. 
Assuming an interaction between the BD scalar field 
and dark matter, one can show that the matter 
dominated era can undergo a transition from a 
decelerated phase to an accelerated phase of 
expansion without using any additional potential 
[17]. It has also been found that the BD scalar field 
alone can cause the acceleration without having to 
consider any quintessence matter or any interaction 
between BD field and dark matter [18]. 


In contrary to the astrophysical observations, many 
of these models predict an ever accelerating 
expansion for the matter dominated universe. To 
account for the recent acceleration, many of the 
models require a very low value of the Brans-Dicke 
parameter w, of the order of unity, but the local 
astronomical experiments require a very high value 
of w [19]. 


In a recent study on f(R,T) gravity, it has been 
shown, that the time evolution of cosmological 
parameters can be correctly predicted without 
having the necessity to incorporate any cosmological 
constant or quintessence field in the model [20]. The 
entire cosmological scenario has been derived from 
the relationf(R,T) =R+aT+T*. Here R and T 
denote the Ricci scalar and trace of the energy- 
momentum tensor respectively. This model shows a 
smooth transition of the universe from a phase of 
decelerated expansion to the phase of accelerated 
expansion, indicated by a change of sign of the 
deceleration parameter (q). Another recent model, 
based on BD theory, has shown this signature flip of 
q, by assuming an inter-conversion between matter 
and dark energy [21]. These models have been 
constructed for a matter dominated universe. 


The purpose of the present work is to study the 
acceleration of the universe in the framework of the 
Generalised Brans-Dicke theory. This version of the 
BD theory was first proposed by Bergmann and 
Wagoner [22, 23]. In the generalised form, the BD 
parameter wis regarded as a function of the BD 
scalar field (@), instead of being a constant. For 
various physical requirements, different functional 
forms of w would originate. Here we have tried to 
find out whether a time-varying w can overcome the 
problems of usual BD models and lead to a smooth 
transition from the phase of decelerated expansion 
to accelerated expansion of the universe. 


In Section 2 of this article, we have written the field 
equations of the generalized BD theory and derived 
the scale factor (a). We have used an ansatz for the 
scalar field (@), expressed in terms of the scale factor 
(a). In Section 3, we have analysed the deceleration 
parameter (q). An expression for w(y)has been 
derived in section 4. In Section 5, we have discussed 
the time evolution of the gravitational constant. 
Section 6 has some discussions regarding the results 
obtained from this model. 
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2. BRANS-DICKE FIELD EQUATIONS AND THE 
SCALE FACTOR 

According to the generalized Brans-Dicke theory, the 
action is given by [24], 


S=f [ + oO 0p" + Ls J—gd*x (1) 


167G 


Here, R denotes the Ricci scalar, Lm denotes the 
matter Lagrangian, g represents the Brans-Dicke 
scalar field andw is a dimensionless coupling 
parameter which is regarded as a function of g here, 
instead of being treated as a constant. 

For a zero-curvature Robertson-Walker space-time, 
the line element is expressed as, 


ds? = dt? —a’(t)f(,6,¢) (2) 
where, f = dr? + r7(d@? + sin?0d¢7) 

Here, a is the scale factor of the universe. Variation of 
action (1) with respect to the metric tensor 


components, given by equation (2), yields the 
following field equations. 


a(e) 242m (ey ate @ 
28+ (0) = sO) -288-2 @ 


In equation (3), p denotes the density of the entire 
matter distribution and 87G = 1. An overhead dot 
stands for differentiation with respect to time. This 
model is based on an assumption that the 
thermodynamic pressure of the cosmic fluid is zero, 
which is quite consistent with the present dust 
universe. To make the above differential equations 
tractable, we have chosen the following relation. 


= a\r _ ay” n 
P = Po (<) = Po (*) = Poa (5) 
0. 

Here adj is the scale factor at the present time. Its 

value has been chosen to be 1, as per convention. 

Considering the conservation of matter, we can 
. _ Pod? _ Po 

write, p= a = a3 

Combining equations (3) and (4) and substituting 

from equation (5), one gets, 


‘i i\2 1 = 
es (n+ 2) (¢) = —_ Po g-(n+3) (6) 


N+2 Po 


The solution of equation (6) gives the following 
expression for the scale factor. 


1 
G= Gem tapo Capel Gen (7) 


2(2+n) Popo 


Here,C, and C, are the integration constants. To 
determine the values of these constants we have 
used the following boundary conditions. 
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1) att =t), a=a) =1 and2)att=ty, d/a=Ho 
Heret, denotes the present epoch. Using these 
boundary conditions, the constants C, and C, are 
found to have the following forms. 


_— y2_ 2Po 
Cy a Ho @o(2+n)(3+n) (8) 
Cz = ano — to (9) 


Using equations (8) and (9) in (7) we get the 
following expression for the scale factor. 


a=[a,+ ain (10) 


= = Po 
f—) a f—) 


Scale factor (a) 
jo) 
a 


0.0 
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ay = 1+ Hy(3 + n)(t o— to) 


apa a, = 2atare-to)? 
, 2 2(2+n)@o 


where, 


The variation of scale factor (a), as a function of 
(t —t)), where t, denotes the present epoch, has 
been shown in Figure 1, for two values of the 
parameter n. 


3. TIME EVOLUTION OF THE DECELERATION 
PARAMETER (q) 

Using equation (10), one obtains the following 
expression for the deceleration parameter. 


4.0x10" 


t- t, (sec) 


Figure 1: Variation of the scale factor (a) as a function of (t — ty), where ty) denotes the present epoch, for 
two values of the parameter n. 
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Figure 2: Variation of the deceleration parameter (q) as a function of (t — t)), where ty denotes the 
present epoch, for two values of the parameter n. 


ai 


a 1 
q= = =a a) (11) 
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ee (2+n)po{PoHo” (3+n)(2+n)-2po} 
2 {(2+n) ~oHo+Po(t-to)}? 


Here, q, =1+7, 


The transition of the universe, from a decelerated to 
an accelerated phase of expansion, indicates clearly 


that the deceleration parameter q was initially 
positive, which became zero at a certain value of t and 
has attained a negative value thereafter, in the recent 
past of the matter dominated era. This implies that the 
present value of the deceleration parameter 
(i. e. the value of gq at t = t)) should be negative. From 
equation (11), we have determined the time (7) at 
which q was zero. It is expressed as, 
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T=t-fath (12) 


where, f, = a. (2 +n) 
0 


and, fo = [eaeeleectatenoenne3 


po?(1+n) 
It is obvious from astrophysical observations 
that T <t ) .This requirement is fulfilled if the 


parameter n satisfies the following condition. 


n<n<ny, (13) 
where, n_ =—2-— [-%, 

PoHo 
and, ne =—2+ [% 

PoHo 


According to BD theory, gis the reciprocal of the 
effective Newtonian’ gravitational constant 4G. 


Therefore, the value of the parameter @) = 
1 


ao 1.498 x 107° N1m?Kg2. The values of other 
constants are, Hy = 2.33x107'8sec1, py = 2.83 x 
10°-77Kg/m and ty = 4.41 x 107’sec, where Hy and 
Po are the values of the Hubble parameter and the 
density of matter (both dark and baryonic) at the 
present epoch (t)) [25]. Using the values of these 
constants in equation (13), the range of values of nis 
obtained as —2.185 <n < —1.814. For the present 
study, we have used two values of n (n = —1.82 and 
—1.88) within this range, chosen in a manner such 
that the time evolutions of the scale factor (a) and the 
deceleration parameter (q) are consistent with the 
observations. The plot of q versus t —t,) has been 
shown in Figure 2 for two values of n. It is evident 
from this plot that the value of q is positive initially, 
becoming zero at a certain value of t—t, and it 
gradually attains a negative value at a value of t for 
which t < t,. It is also observed in the same figure 
that the value of the deceleration parameter (q) 
remains negative even at the current epoch and 
thereafter. One may choose different values ofn, 
within the obtained range, to just fine tune the value 
of t—t, at which the signature flip of q occurs, 
without affecting the nature of the plot. 


4. FUNCTIONAL FORM OF THE BRANS-DICKE 
PARAMETER w(@) 

Combining equations (3), (4), (5) and the expression 
of p = p,/a?, one obtains the following expression for 
the dimensionless parameter w(@) of BD theory. 


2 


w(y) = 3 fac +n) - 20-6 (2) "| (14) 


n2 a 


Using equation (10) in (14) we get, 


w(~) = Slo, - a2] (15) 
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where, @, = 6(1+7n) 
2po(2+n)(3+n)(9/~o)*3/” 


~ Ho? @o0(2+N)(3+N)+2po{-1+(~/Po)1t3/} 


and, W 


Equation (15) shows that for n = —1.85, we have 
w ~ —1.49 for an extremely large value of ~ and, for 
an extremely small value of ~, w ~ —1.54. Since n is 
negative, @ must be decreasing with time, in 
accordance with equation (5). The above analysis 
therefore shows that w becomes more and more 
negative with time. 


5. TIME EVOLUTION OF THE GRAVITATIONAL 
CONSTANT (G) 

Since mis found to be decreasing with time, the 
gravitational constant G(= 1/g) must be increasing 
with time. This result is consistent with other recent 
studies [26-28]. Combining equations (5) and (10), we 
get the following expression for G. 


—n 


1 1 1 exis 
G=—=—a™" =—[a, + a,|3" 16 

a a = 1+ a)] (16) 
where, a, = 1+4H)(3+n)(t — to) 


and, a, = 2oetmte=to)” 
: 2 2(2+N)Po 


6. CONCLUSIONS 

The present article explains the accelerated 
expansion of the universe in the framework of 
generalized Brans-Dicke theory. It clearly 
demonstrates that a time-dependent w can indeed 
generate a signature flip in the deceleration 
parameter (q) in a matter dominated universe, unlike 
many other models based on Brans-Dicke theory. The 
plot, that shows the time dependence of gq, clearly 
indicates a smooth change of the deceleration 
parameter from a positive to a negative value, in the 
recent past. Therefore, it is proved that the 
generalized BD theory is capable of explaining the 
observed cosmic acceleration, without taking into 
consideration any component that represents dark 
energy. This theory has an advantage due to the fact 
that the scalar field is already there as an integral part 
of the theory. Here, the dimensionless coupling 
parameter (w) is found to have a negative value. Its 
functional dependence upon the scalar field (~) has 
played an important role in making the effective 
pressure negative. 

However, this theoretical study has its own 
shortcomings. Although one gets the correct nature of 
time evolution of the cosmological parameters, this 
model has failed to generate the correct value of the 
deceleration parameter for t = ty. For better results, 
one may try different empirical expressions for 9, 
instead of the one represented by equation (5). Apart 
from that, the BD parameter (w) has a low negative 
value and this is in contradiction to the local 
astronomical requirement of a high value of w. This 
article, however, provides a direction in which one 
may carry out theoretical studies so that a suitable 
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form of w(g) can be obtained, which would explain 
the cosmic acceleration properly and, at the same 


time, 


generate a value of w that is consistent with the 


local astronomical observations. 
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